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Abstract 

Overlap-free words are words over the binary alphabet A = {a, b} that do not 
contain factors of the form xvxvx, where x G A and v £ A* . We analyze the 
asymptotic growth of the number u n of overlap-free words of length n as n — > oo. We 
obtain explicit formulas for the minimal and maximal rates of growth of u n in terms 
of spectral characteristics (the lower spectral radius and the joint spectral radius) 
of certain sets of matrices of dimension 20 x 20. Using these descriptions we provide 
new estimates of the rates of growth that are within 0.4% and 0.03% of their exact 
values. The best previously known bounds were within 11% and 3% respectively. We 
then prove that the value of u n actually has the same rate of growth for "almost all" 
natural numbers n. This "average" growth is distinct from the maximal and minimal 
rates and can also be expressed in terms of a spectral quantity (the Lyapunov 
exponent). We use this expression to estimate it. In order to obtain our estimates, 
we introduce new algorithms to compute spectral characteristics of sets of matrices. 
These algorithms can be used in other contexts and are of independent interest. 

Key words: Overlap-free words, Combinatorics on words, Joint spectral radius, 
Lyapunov exponent. 



Preprint submitted to Elsevier 



1 Introduction 



Binary overlap-free words have been studied for more than a century. These 
are words over the binary alphabet A = {a, b} that do not contain factors 
of the form xvxvx, where x G A and v G A*. For instance, the word baabaa 
is overlap free, but the word baabaab is not, since it can be written xuxux 
with x = b and u = aa. See [1] for a recent survey. Thue [22, 23] proved in 
1906 that there are infinitely many overlap- free words. Indeed, the well-known 
Thue- Morse sequenc^ZI is overlap-free, and so the set of its factors provides an 
infinite number of different overlap-free words. The asymptotics of the number 
u n of such words of a given length n was analyzed in a number of subsequent 
contribution^]. The number of factors of length n in the Thue-Morse sequence 
is proved in [6] to be larger than 3n, thus providing a linear lower bound on 

u n > 3n. 

The next improvement was obtained by Restivo and Salemi [20]. By using 
a certain decomposition result, they showed that the number of overlap-free 
words grows at most polynomially: 

u n < Cn r , 

where r = log(15) ~ 3.906. This bound has been sharpened successively by 
Kfoury [11], Kobayashi [12], and finally by Lepisto [13] to the value r = 
1.37. One could then suspect that the sequence u n grows linearly. However, 
Kobayashi [12] proved that this is not the case. By enumerating the subset of 
overlap-free words of length n that can be infinitely extended to the right he 
showed that u n > C n 1 - 155 and so we have 

dn 1 ' 155 < u n < C 2 n L37 . 

It is worth noting that the sequence u n is 2-regular, as shown by Carpi [7]. 
On Figure 1(a) we show the values of the sequence u n for 1 < n < 200 and on 
Figure 1(b) we show the behavior of logu n /logn for larger values of n. One 
can see that the sequence u n is not monotonic, but is globally increasing with 
n. Moreover the sequence does not appear to have a polynomial growth since 
the value \ogu n / \ogn does not seem to converge. In view of this, a natural 
question arises: is the sequence u n asymptotically equivalent to n r for some r ? 
Cassaigne proved in [8] that the answer is negative. He introduced the lower 
and the upper exponents of growth: 



1 The Thue-Morse sequence is the infinite word obtained as the limit of 6 n (a) for 
n — > oo with 9(a) = ab, 0(b) = ba; see [8]. 

2 The number of overlap-free words of length n is referenced in the On-Line Ency- 
clopedia of Integer Sequences under the code A007777; see [21]. The sequence starts 
1, 2, 4, 6, 10, 14, 20, 24, 30, 36, 44, 48, 60, 60, 62, 72,... 
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(a) (b) 

Fig. 1. The values of u n for 1 < n < 200 (a) and logtt n /logn for 1 < n < 10000 
(b). 

a = supjr BC > 0, u n > CV}, (1) 
ft = mi[r \3C > 0,u n < Cn r }, 

and showed that a < ft. Cassaigne made a real breakthrough in the study 
of overlap-free words by characterizing in a constructive way the whole set of 
overlap-free words. By improving the decomposition theorem of Restivo and 
Salemi he showed that the numbers u n can be computed as sums of variables 
that are obtained by certain linear recurrence relations. These relations are 
explicitly given in the next section and all numerical values can be found in 
Appendix A. As a result of this description, the number of overlap-free words 
of length n can be computed in logarithmic time. For the exponents of growth 
Cassaigne has also obtained the following bounds: a < 1.276 and ft > 1.332. 
Thus, combining this with the earlier results described above, one has the 
following inequalities: 

1.155 < a < 1.276 and 1.332 < ft < 1.37. (2) 

In this paper we develop a linear algebraic approach to study the asymptotic 
behavior of the number of overlap-free words of length n. Using the results 
of Cassaigne we show in Theorem 2 that u n is asymptotically equivalent to 
the norm of a long product of two particular matrices Aq and A\ of dimension 
20 x 20. This product corresponds to the binary expansion of the number n — 1. 
Using this result we express the values of a and ft by means of certain joint 
spectral characteristics of these matrices. We prove that a = log 2 p(A ,Ai) 
and ft = log 2 p(A , Ai), where p and p denote, respectively, the lower spectral 
radius and the joint spectral radius of the matrices Aq,A± (we define these 
notions in the next section). In Section 3, we estimate these values and we 
obtain the following improved bounds for a and ft: 

1.2690 < a < 1.2736 and 1.3322 < ft < 1.3326. (3) 
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Our estimates are, respectively, within 0.4% and 0.03% of the exact values. 
In addition, we show in Theorem 3 that the smallest and the largest rates of 
growth of Uyi £1X6 effectively attained, and there exist positive constants Ci, C<i 
such that C\ n a < u n < C 2 n 13 for all n G N. 

Although the sequence u n does not exhibit an asymptotic polynomial growth, 
we then show in Theorem 5 that for "almost all" values of n the rate of growth 
is actually equal to a = log 2 p(A ,Ai), where p is the Lyapunov exponent 
of the matrices. For almost all values of n the number of overlap-free words 
does not grow as n a , nor as n 13 , but in an intermediary way, as n° ' . This means 
in particular that the value l °^ u ^ converges to a as n — > oo along a subset of 
density 1. We obtain the following bounds for the limit a, which provides an 
estimation within 0.8% of the exact value: 



1.3005 < a < 1.3098. 



These bounds clearly show that a < a < (3. 



To compute the exponents a and a we introduce new efficient algorithms for 
estimating the lower spectral radius p and the Lyapunov exponent p of matri- 
ces. These algorithms are both of independent interest as they can be applied 
to arbitrary matrices. 



Our linear algebraic approach not only allows us to improve the estimates of 
the asymptotics of the number of overlap-free words, but also clarifies some 
aspects of the nature of these words. For instance, we show that the "non 
purely overlap-free words" used in [8] to compute Ufi £1X6 asymptotically neg- 
ligible when considering the total number of overlap-free words. 



The paper is organized as follows. In the next section we formulate and prove 
the main theorems (except for Theorem 2, whose proof is quite technical and is 
given in Appendix B). Then in Section 3 we present algorithms for estimating 
the joint spectral radius, the lower spectral radius, and the Lyapunov exponent 
of linear operators. Applying them to those special matrices we obtain the 
estimations for a, (3 and a. In the appendices we write explicit forms of the 
matrices and initial vectors used to compute u n , we give a proof of Theorem 2 
and present the results of our numerical algorithms. 
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2 The asymptotics of the overlap-free words 



In the sequel we use the following notation: M. d is the <i-dimensional space, 
inequalities x > and A > mean that all the entries of the vector x (respec- 
tively, of the matrix A) are nonneg ative. We denote M d + = {i6 R d , x > 0}, by 

| x | we denote a norm of the vector x G M, d , and by || • || any matrix norm. In par- 

d 

ticular, \x\i = \ x i\i ll^lli = su Pi;ch=i \Ax\ = maXj=i i ... t iX)i=i We write 

i=l 

1 for the vector (1, . . . , 1) T G IR d , p{A) for the spectral radius of the matrix A, 
that is, the largest magnitude of its eigenvalues. If A > 0, then there is a vec- 
tor v > such that Aw = p(A)u (the so-called Perron- Frobenius eigenvector). 
For two functions fi, f% from a set y to IR + the relation fi(y) x /^(y) means 
that there are positive constants Ci,C 2 such that Cif\{y) < f%(y) < C 2 fi(y) 
for all j/6V. 

To compute the number w n of overlap-free words of length n we use several 
results from [8] that we summarize in the following theorem: 



Theorem 1 Let F ,Fi G M 30x30 , and let w,y 3 , - - . , 2/15 G M 30 &e as given in 
Appendix A. For n > 16, let y n be the solution of the following recurrence 
equations 

V2n = F y n 

(4) 

y 2n +i = F x y n . 

Then, for any n > 16, the number of overlap-free words of length n is equal to 
w T y n ^i. 

It follows from this result that the number u n of overlap-free words of length 
n > 16 can be obtained by first computing the binary expansion ■ ■ • d\ of 
n — 1, i.e., n — 1 = X^=o ^j+i^ '\ an d then defining 

u n = w T F dl ■ ■ ■ F dk 4 y m (5) 

where m = d^-s + dk-2^ + <ifc-i2 2 + d^ 3 . To arrive at the results summarized 
in Theorem 1, Cassaigne builds a system of recurrence equations allowing the 
computation of a vector U n whose entries are the number of overlap-free words 
of certain types (there are 16 different types). These recurrence equations also 
involve the recursive computation of a vector V n that counts other words of 
length n, the so-called "single overlaps". The single overlap words are not 
overlap-free, but have to be computed, as they generate overlap-free words of 
larger lengths. We now present the main result of this section which improves 
the above theorem in two directions. First we reduce the dimension of the 
matrices from 30 to 20, and second we prove that u n is given asymptotically 
by the norm of a matrix product. The reduction of the dimension to 20 has 
a straightforward interpretation: when computing the asymptotic growth of 
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the number of overlap-free words, one can neglect the number of "single over- 
laps" V n defined by Cassaigne. We call the remaining words purely overlap-free 
words, as they can be entirely decomposed in a sequence of overlap-free words 
via Cassaigne's decomposition (see [8] for more details). 

Theorem 2 Let Aq, A\ G M^ 0x20 be the matrices defined in Appendix A (Equa- 
tion A. 3), let || • || be a matrix norm, and let A(n) : N — > IR+ 0x20 be defined as 
A{n) = A dl ■ ■ ■ Ad k with d^. . ,d\ the binary expansion of n — 1. Then, 

u n ^\\A{n)\\. (6) 

Observe that the matrices F , Fi in Theorem 1 are both nonnegative and hence 
possess a common invariant cone K = R+. We say that a cone K is invariant 
for a linear operator B if BK C K. All cones are assumed to be solid, convex, 
closed, and pointed. We start with the following simple result proved in [18]. 

Lemma 1 For any cone K C M d ; for any norm \ ■ \ in M. d and any matrix 
norm || • || there is a homogeneous continuous function 7 : K — > M + positive on 
mtK such that for any x G inti^ and for any matrix B that leaves K invariant 
one has 

7(x)||i?|| • \x\ < \Bx\ < \\B\\ ■ \x\. 

70) 

Corollary 1 Let two matrices Aq,A± possess an invariant cone K C M. d . 
Then for any x G inti^ we have ■ ■ ■ Ad k x\ x \\Afa ■ ■ -A^W for all k and 
for all indices di, . . . , d k G {0, 1}. 

In view of Corollary 1 and of Eq. (5), Theorem 2 may seem obvious, at least if 
we consider the matrices Fi instead of A%. One can however not directly apply 
Lemma 1 and Corollary 1 to the matrices Aq, A\ or to the matrices Fq, Fi 
because the vector corresponding to x is not in the interior of the positive 
orthant, which is an invariant cone of these matrices. To prove Theorem 2 
we construct a wider invariant cone of A and A\ by using special properties 
of these matrices. We detail the construction of this cone in the proof given 
in Appendix B. Theorem 2 allows us to express the rates of growth of the 
sequence u n in terms of norms of products of the matrices Aq, A\ and then 
to use joint spectral characteristics of these matrices to estimate the rates of 
growth. More explicitly, Theorem 2 yields the following corollary: 

Corollary 2 Let Aq, A\ G M 2 h 0x20 be the matrices defined in Appendix A and 
let A(n) : N —> Rf x20 be defined as A(n) = A dl ■ ■ ■ A dk with d k . . . d 1 the 
binary expansion of n — 1 . Then 

l ^^-\og\\A(n)\\ 1 / k -> as n — 00. (7) 
log n 
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Proof. Since log fc 2 11 — > 1 as n — > oo, we have 



lim 



n— »oo \ lu 62 



log 2 u n _ log 2 ll^ dl ---A dfc | 
log 2 n fc 



i lQg 2 U„-log 2 ||A dl ...A d J| = ^ log^Un-ll^-^H-l) 



n— >oo 



By Theorem 2 the value \og 2 (u n ■ \\A dl ■ ■ ■ A dk is bounded uniformly over 
n G N, hence it tends to zero, being divided by k. □ 
We first analyze the smallest and the largest exponents of growth a and f3 
defined in Eq. (1). For a given set of matrices E = {A\, . . . , A m } we denote 
by p and p its lower spectral radius and its joint spectral radius: 



p(£)=lim mm \\A dl ■ ■ ■ A dk f/\ (8) 

k— >oo di,...,d k £{l,...,m} 

/3(E) = lim max • ■ ■ A d J 1/fc . 

fe— >oo ai,...,Ofc6{l,...,m} 

Both limits are well-defined and do not depend on the chosen norm. Moreover, 
for any product A dl ■ ■ ■ A dk we have 

P<p(A dl .-.A dk y/ k <p (9) 

(see [16,5] for surveys on these notions). 

Theorem 3 For k > 1, let a k = min and 3 k = max ^ SJ±IL - 
Then 

a = lim a k = log 2 p(A , A) and j3 = lim /? fc = log 2 p(A , A x ), (10) 

where the matrices Aq,A x are defined in Appendix A. Moreover, there are 
positive constants G\,Gi such that 

C\ < min u n n~ a and C\ < max u n n~^ < Ci (11) 

2 k ~ 1 <n<2 k 2 k ~ 1 <n<2 k 

for all k G N. 

The proof of this theorem is based on the following auxiliary result taken 
from [18]. For a given set of indices {ii, . . . , i p } C {1, . . . , d}, p > 1 we call 
the subspace L^,...,^ = {x G M d , x^ — ■ ■ ■ — Xi — 0} a coordinate plane. 

Proposition 1 [18] Let A$,Ai be matrices with a common invariant cone. 
Then there is a positive constant c\ such that 

max \\A dl ■ • • A dk \\ > cip k and min \\A dl ■ ■ ■ A d J > c 1 p k , k G N. 

di,...,d k di,...,d k 
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If, moreover, these matrices have no common invariant subspace among the 
coordinate planes, then there is a positive constant c 2 such that 

max \\A dl ■ ■ ■ A dk \\ < c 2 p k , k e N. 

di,...,dk 



Proof of Theorem 3. The equalities in Eq. (10) follow immediately from Corol- 
lary 2 and the definitions given on Eq. (8). To prove the inequalities given at 
Eq. (11) we apply Proposition 1 for our matrices Aq, A\ that have an invariant 
cone Theorem 2 yields 

u n n~ a x \\A dl ■ ■ -A dk \\2- ak = \\A dl ■ ■ ■ A dk \\p~ k . 

Taking the minimum over n = 2 fc_1 + 1, . . . , 2 k and invoking Proposition 1 we 
conclude that min u n n~ a > C\. The same with the inequality max u n n~P > 

2 k - 1 <n<2 k 2 k - 1 <n<2 k 

C\. To prove the upper bound in Eq. (11) we note that the matrices A ,Ai 
have no common invariant subspaces among the coordinate planes (to see this 
observe, for instance, that (A Q + Ai) 5 has no zero entry). □ 



Corollary 3 There are positive constants C\,C 2 such that 



C t n a < u n < C 2 n p , n e N. 



In the next section we will see that a < f3. In particular, the sequence u n does 
not have a constant rate of growth, and the value l ° gUn does not converge as 

' log n 

n — ► 00. This was already noted by Cassaigne in [8]. Nevertheless, it appears 
that the value actually has a limit as n — > 00, not along all the natural 
numbers n G N, but along a subsequence of N of density 1. In other terms, the 
sequence converges with probability 1. The limit, which differs from both a 
and P can be expressed by the so-called Lyapunov exponent p of the matrices 
Aq, A\. To show this we apply the following result proved by Oseledets in 1968. 
For the sake of simplicity we formulate it for two matrices, although it can be 
easily generalized to any finite set of matrices. 



Theorem 4 [14] Let A Q , A\ be arbitrary matrices and d±, c^, • • ■ be a sequence 
of independent random variables that take values and 1 with equal probabil- 
ities 1/2. Then the value \\A dl • • ■ A dk \\ 1 ' k converges to some number p with 



probability 1. This means that for any e > we have P 
s) — > as k — > 00. 



P 



> 



8 



The limit p in Theorem 4 is called the Lyapunov exponent of the set {A Q , A{\. 
This value is given by the following formula: 

p(A ,A 1 )= lim ( J] \\A dl ---M\\ 1/k ) 1/2 (12) 

5 7 



(for the proof see, for instance, [19]). To understand what this gives for the 
asymptotics of our sequence u n we introduce some further notation. Let V be 
some property of natural numbers. For a given k G N we denote 

P k (V)= 2~ (fc - 1) Card{n e {2 fc_1 + 1, . . . , 2*}, n satisfies V }. 

Thus, Pfe is the probability that the integer n uniformly distributed on the 
set {2 + 1, . . . , 2 k } satisfies V. Combining Proposition 2 and Theorem 4 we 
obtain 

Theorem 5 There is a number a such that for any e > we have 

logu n 



Pk 



log 72 



a 



> e — > as k — ► oo. 



Moreover, a = log 2 p, where p is the Lyapunov exponent of the matrices 
{Ao,Ai} defined in Appendix A. 

Thus, for almost all numbers n G N the number of overlap-free words u n has 
the same exponent of growth a = log 2 p. If positive a and b are large enough 
and a < b, then for a number n taken randomly from the segment [a, b] the 
value : ° g is close to a. Let us recall that a subset A C N is said to have 

log n 

density 1 if ^ Card jr < n, r G «4.| — > 1 as n — > oo. We say that a sequence 
/„ converges to a number / along a set of density 1 if there is a set A C N of 
density 1 such that lim /„ = /. Theorem 5 yields 

n— >oo,n€_4 



Corollary 4 TTie value l ° gUn converges to a along a set of density 1. 

17 log n ^ z j & 

Proof. By Theorem 5 there exists ki G M such that at least half of the natural 
numbers n < 2 hl satisfy the inequality — cr < |. Denote the set of such 

numbers n < 2 kl by A\. Further, there exists k 2 > k\ such that at least ~ 



of the numbers n G {2 A 



2 fc2 } satisfy the inequality 



log u n 
log n 



cr 



< 



i 



Denote the set of such numbers by A 2 - The remainder is by induction: having 



a number k r . 
numbers n G 



we take a number k r > k r -\ such that at least 



-i 



{2 h 



+ 1, 



2 fer } satisfy the inequality 



log Un 

log n 



(7 



2' 

< 



of the 
and 



denote the set of such numbers by A T . The set A = U rS N^4r has density 1 and 
l ° gUn tends to a along A as n — >• oo. □ 
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3 Estimations of the exponents 



Theorems 2 and 5 reduce the problem of estimating the exponents of growth 
of u n to computing joint spectral characteristics of the matrices A and A\. 
In order to estimate the joint spectral radius we use a modified version of 
the "ellipsoidal norm algorithm" [3]. For the lower spectral radius and for the 
Lyapunov exponent we present new algorithms, which seem to be relatively 
efficient for nonnegative matrices. The results we obtain can be summarized 
in the following theorem: 



Theorem 6 

1.2690 < a < 1.2736 

1.3322 < (3 < 1.3326 (13) 
1.3005 < a < 1.3098 

In this section we also make (and give arguments for) the following conjecture: 
Conjecture 1 



= log 2 ^p(A Q A 1 ) = 1.3322.... 
3. 1 Estimation of (3 and the joint spectral radius 



By Theorem 3 to estimate the exponent (3 one needs to estimate the joint 
spectral radius of the set {v4o,Ai}. A lower bound for p can be obtained by 
applying inequality (9). Taking k — 2 and d\ — 0, = 1 we get 

p > [pOMi)] = 2.5179, (14) 

and so (3 > log 2 2.5179 > 1.3322 (this lower bound was already found in [8]). 
Upper bounds for the joint spectral radius of sets of matrices E = {A±, . . . , A m } 
are usually derived from the following simple inequality 

p < max \\A dl ---A dk \\ l '\ (15) 

di,...,d k e {1,... ,m} 

which holds for every k > 1 and converges to p as k — > oo. This, at least 
theoretically, gives arbitrarily sharp estimations for p. However, in our case, 
due to the size of the matrices A ,Ai, this method leads to computations 
that are too expensive even for relatively small values of k. Faster convergence 
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can be achieved by finding an appropriate norm. To do this we use the so- 
called ellipsoidal norm: \\A\ \p = max x \j xT ^Ax ^-- wnere -P is a positive definite 
matrix. This is the matrix norm induced by the vector norm \x\p — (x T Px) 1 ^ 2 . 
The crucial idea is that the optimal P, for which the right hand side in (15) for 
k — 1 is minimal, can be found by solving a simple semidefinite programming 
problem. This algorithm can be iterated using the relation = In 

the sequel we denote £ fc = { A dl ■ ■ ■ A dk , 1 < di < m,i = 1, . . . , k j. Thus 

one can consider the set E fc as a new set of matrices, and approximate its joint 
spectral radius with the best possible ellipsoidal norm. In Appendix C we give 
an ellipsoidal norm such that each matrix in S 14 has a norm smaller than 
2.5186 14 . This implies that p < 2.5186, which gives (3 < 1.3326. Combining 
this with the inequality j3 > 1.3322 we complete the proof of the bounds for 
/3 in Theorem 6. 

We have not been able to improve the lower bound of Eq. (14). However, 
the upper bound we obtain is very close to this lower bound, and the upper 
bounds obtained with an ellipsoidal norm for E fc get closer and closer to this 
value when k increases. Moreover, it has already been observed that for many 
sets of matrices for which the joint spectral radius is known exactly, and in 
particular matrices with nonnegative integer entries, there always is a product 
that achieves the joint spectral radius, i.e., a product A e S* such that p = 
p(AY 1 ^ [9]. For these reasons, we conjecture that the exponent j3 is actually 
equal to the lower bound. 



3.2 Estimation of a and the lower spectral radius 



An upper bound for p(A Q ,Ai) can be obtained using Eq.(9) for k = 1 and 
d\ = 0. We have 

a = log 2 (p) < log 2 (A ) = 1.276... (16) 

This bound for a was first derived in [8]. This bound is however not optimal. 
Taking the product A}°y4 (i.e., k = 11 in inequality (9)), we get a better 
estimate: 

a < log 2 [(p^fAo) 1 / 11 ] = 1.2735... (17) 

One can verify numerically that this product gives the best possible upper 
bound among all the matrix products of length k < 14. 

We now estimate a from below. The problem of approximating the lower spec- 
tral radius is NP-hard [4] and to the best of our knowledge, no algorithm is 
known to compute p. Here we propose two new algorithms. We first consider 
nonnegative matrices. As we observed above, for any k we have p(S fc ) = p fc (S). 
Without loss of generality it can be assumed that the matrices of the set E do 
not have a common zero column. Otherwise, by suppressing this column and 
the corresponding row we obtain a set of matrices of smaller dimension with 
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the same lower spectral radius. The vector of ones is denoted by 1. 

Theorem 7 Let E be a set of nonnegative matrices that do not have any 
common zero column. If for some r G < t G N, there exists x G M. d 

satisfying the following system of linear inequalities 

B(Ax - rx) > 0, WB G E s , A G £* 

(18) 

x > 0, (x, 1) = 1, 

f/ien p(E) > r 1 /*. 

Proof. Let x be a solution of (18). Let us consider a product of matrices 
A k ...Ax G E fe * : Aj G E*. We show by induction in k that A k . . . A x x > 
r k ~ 1 A k x : For k = 2, we have A2(A x x — rx) = CB(A x x — rx) > 0, with 
£ G E S ,C G E*~ s . For k > 2 we have A fe ...Aix = A fc v4 fc _i . . . A x x > 
r k ~ 2 AkAk- X x > r k ~ 1 AkX. In the last inequality the case for k = 2 was 
reused. Hence, \\A k . . . A x \\ = l T A k . . . A X 1 > r^ 1 ! 7 A k x > r k C, where C = 
(minfe l T A k x) jr > 0. The last inequality holds because A k x = 0, together 
with the first inequality in (18), imply that —rBx = for all B G E s , which 
implies that all B G E s have a common zero column. This is in contradiction 
with our assumption because the matrices in E s share a common zero column 
if and only if the matrices in E do. □ 
We were able to find a solution to the linear programming problem (18) 
with r = 2.41, t = 16, s — 6. Hence we get the following lower bound: 
a > logr > 1.2690. The corresponding vector x is given in Appendix D. 
This completes the proof of Theorem 6. 

Theorem 7 handles nonnegative matrices, and we propose now a way to gener- 
alize this result to arbitrary real matrices. The idea is to lift the matrices to a 
larger vector space, so that all the matrices share an invariant cone. This kind 
of lifting is rather classical and is known under several names in the literature 
as for instance semidefinite lifting or symmetric algebras [2,17,15]. The idea 
is to consider the matrices A^ G E as linear operators acting on the cone of 
positive semidefinite matrices S as S — > AJSA4. It is not difficult to prove 
that the lower spectral radius of this new set of linear operators is equal to 
p(E) 2 . We use the notation A y B to denote that the matrix A — B is positive 
semidefinite. Recall that A y <^ Vy, y 7 Ay > 0. 

Theorem 8 Let E be a set of matrices in M. dxd and s < t G N. Suppose that 
there are r > and a symmetric matrix S y such that 

B*(A*SA - rS)B y \/A G E*, B G E s 

(19) 

SyO 
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then p(S) > r l l 2t . 

Proof. The proof is formally similar to the previous one: Let S be a solution 
of (19). We denote My. the product A\ . . . Af.. It is easy to show by induction 
that MlSM k y r k ~\AlSA k ). This is obvious for k = 2 for similar reasons as 
in the previous theorem, and for k > 2, if, by induction, 

Vy, y*M* k _ 1 SM k _ 1 y > r k ~ 2 y* A^SA^y, 

then, with y = A k x, for all x, 

x*M* k SM k x > r k - 2 x*A* k A* k _ 1 SA k _ 1 A k x > r fc_1 z* A* k SA k x. 

Thus, 

f x*M k *SM k x \ y ^ fc _ lgu f x*A* k SA k x \ 
I x*Sx ) ~ I J 

Finally, ||M fe || 5 > r k ' 2 C, where C is a constant. □ 
For a given r > the existence of a solution 5 can be established by solving 
the semidefmite programming problem (19), and the optimal r can be found 
by bisection in logarithmic time. 

3.3 Estimation of a and the Lyapunov exponent 

The exponent of the average growth a is obviously between a and (3, so 
1.2690 < a < 1.3326. To get better bounds we need to estimate the Lya- 
punov exponent p of the matrices Aq, A\. The first upper bound can be given 
by the so-called 1-radius p\. 

Pl = lim(2- fc £ \\A dl ---A dk \\) 1/k . 

fe ^°° v d lt ...,d k J 

For matrices with a common invariant cone we have p\ = \p(A$ + Ai) [18]. 
Therefore, in our case p\ = \p(Aq + Ai) = 2.479.... This exponent was first 
computed in [8], where it was shown that the value YJjZo u j is equivalent 
to n 71 , where 77 = 1 + log 2 pi = 2.310.... It follows immediately from the in- 
equality between the arithmetic mean and the geometric mean that p < p\. 
Thus, a < r\. In fact, as we show below, a is strictly smaller than 77. We are 
not aware of any approximation algorithm for the Lyapunov exponent, except 

by application of Definition (12). It is easily seen that for any k the value 

1 

r k = ( II 1 1 Aii ' ' ' Ad k \\ J gives an upper bound for p, that is p < r k for 

any fceN. Since r k — > p as k — > 00, we see that this estimation can be arbitrar- 
ily sharp for large k. But for the dimension 20 this leads to extensive numerical 
computations. For example, for the norm || • ||i we have r2o = 2.4865, which is 
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even larger than p\. In order to obtain a better bound for p we state the follow- 
ing results. For any k and x £ M. d we denote Pk(x) — ( Yl \A dl ■ ■ ■ A dk x\ 

di,...,d k 

and m k = sup Pk{x)- 

x>0,|x|=l 

Proposition 2 LetA ,A 1 be nonnegative matrices in ~R d . Then for any norm 
| • | and for any k > 1 we have p < (m^) 1 ^. 

Proof. For any x > and for any n, k £ N we have pfc +n (x) < m^ n {x). 

r I 1 /tfc 

Therefore, Ptk( x ) — (^fc)*- Whence, lim ptk(x) < f 771 !) 1 — (^fc) 1//fc - On 
the other hand, by Corollary 1 for x > we have r n x [^(a;)] 1 /™, and conse- 

r -il/tfc . 

quently lim ptk{x) — » p as t — »■ oo. Thus, p < (m^) ' . □ 

i— >oo L 

Proposition 3 Lei Ao, Ai fre nonnegative matrices in ~R d i/iai do not /iai>e 

common invariant subspaces among the coordinate planes. If p < p, then p < 
Pi- 



Proof. Let v* be the eigenvector of the matrix |L4o-|-A*j corresponding to its 
largest eigenvalue p\. Since the matrices have no common invariant coordinate 
planes, it follows that > 0. Consider the norm |x| = (x, v*) on R^. Take 
some fc > 1 and y £ Ri, |y| = (y, u*) = 1, such that Pfc(y) = mj.. We have 

m k = p k (y) < 2~ k \M ■ ■ ■ A dk v\ = 2 ~ k J2 { A di ' • • A d h V, v.) 

di,...,d k d 1 ,...,d k 



Thus, m fc < p^, and the equality is possible only if all 2 fc values \A dl ■ ■ ■ A dk y\ 
are equal. Since p < p, there must be a such that the inequality is strict. 
Thus, m k < p\ for some k, and by Proposition 2 we have p < (mk) 1 ' k < Pi-Q 



We are now able to estimate p for the matrices Aq,Ai. For the norm \x\ = 
(x,v*) used in the proof of Proposition 3 the value — -|log 2 m fc can be found 
as the solution of the following convex minimization problem with linear con- 
straints: 

di,...,d k e{o,i} \ / (20) 

s.t. x > 0, (x, v*) = 1. 

The optimal value of this optimization problem is equal to — \ log 2 nik, which 
gives un upper bound for a = log 2 p (Proposition 2). Solving this problem for 
k = 12 we obtain a < 1.3098. We finally provide a theorem that allows us to 
derive a lower bound on a. The idea is identical to the one used in Theorem 7, 



11 



but transposed to the Lyapunov exponent. 



Theorem 9 Let E be a set of nonnegative matrices that do not have any 
common zero column. If for some rj G < t G N, there exists x G 

satisfying the following system of linear inequalities 

B(A iX - nx) > 0, \/B G E s , Ai G E*, 
x > 0, (x, 1) = 1, 

toen p(E) > n^ 1 ^ 2 '). 

The proof is similar to the proof of Theorem 7 and is left to the reader. Also, a 
similar theorem can be stated for general matrices (with negative entries), but 
involving linear matrix inequalities. Due to the number of different variables 
r*j, one cannot hope to find the optimal x with SDP and bisection techniques. 
However, by using the vector x computed for approximating the lower spectral 
radius (given in Appendix D), with the values s — 8, t — 16 for the parameters, 
one gets a good lower bound for a : o > 1.3005. 



4 Conclusions 



The goal of this paper is to precisely characterize the asymptotic rate of growth 
of the number of overlap-free words. Based on Cassaigne's description of these 
words with products of matrices, we first prove that these matrices can be 
simplified, by decreasing the state space dimension from 30 to 20. This im- 
provement is not only useful for numerical computations, but allows to char- 
acterize the overlap-free words that "count" for the asymptotics: we call these 
words purely overlap free, as they can be expressed iteratively as the image of 
shorter purely overlap free words. 

We have then proved that the lower and upper exponents a and (5 defined by 
Cassaigne are effectively reached for an infinite number of lengths, and we have 
characterized them respectively as the logarithms of the lower spectral radius 
and the joint spectral radius of the simplified matrices that we constructed. 
This characterization, combined with new algorithms that we propose to ap- 
proximate the lower spectral radius, allow us to compute them within 0.4%. 
The algorithms we propose can of course be used to reach any degree of accu- 
racy for j3 (this seems also to be the case for a and a, but no theoretical result 
is known for the approximation of the lower spectral radius). The computa- 
tional results we report in this paper have all been obtained in a few minutes 
of computation time on a standard PC desktop and can therefore easily be 
improved. 
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Finally we have shown that for almost all values of n, the number of overlap- 
free words of length n do not grow as n a , nor as n", but in an intermediary 
way as n a , and we have provided sharp bounds for this value of a. 
This work opens obvious questions: Can joint spectral characteristics be used 
to describe the rate of growth of other languages, such as for instance the 
more general repetition free languages ? The generalization does not seem to 
be straightforward for several reasons: first, the somewhat technical proofs of 
the links between u n and the norm of a corresponding matrix product take 
into account the very structure of these particular matrices, and second, it 
is known that a bifurcation occurs for the growth of repetition-free words: 
for some members of this class of languages the growth is polynomial, as for 
overlap-free words, but for some others the growth is exponential [10], and 
one could wonder how the joint spectral characteristics developed in this pa- 
per could represent both kinds of growth. 
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A Numerical values 



We introduce the following auxiliary matrices. For the sake of simplicity our 
notation do not follow exactly those of [8]. 



(o 











1 


2 






f° 








12 1 


\ 




( 




















2 4 2^ 





1 


1 1 


1 

































1 


1 





1 1 

















1 


1 













1 


10 



















1 1 


1 1 




























1 


1 







n 

u 





1 


1 











1 


2 


1 


o 


o 


o 


o 




o 













n 
u 


o 


o 


() 


o 
























,B X = 














,Ci = 
















o 


1 


1 


o 


o 


() 


o 




o 
















1 


o 


o 


1 













































1 


1 





























1 
















1 






















2 































1 
















1 































{ V 




v 








c 


1 




[° 

















o y 








(o 











o) 




/ 

(o 





0^ 





























































































1 


1 

















f 1 1 1 1 


) 

























1 


1 





































































1 1 






















































































,c 2 = 










G\ = 


110 


































































































1 




































































































1 


J 













































2 
































1 



















1 

































1° 















v° 





o o) 

























Now, defining 





OlOxlO 


c 2 


010x5^ 




( Dl 


Bi 


0l0x5 


bA 


Di 


By 


0l0x5 


B 2 




OlOxlO 


Ci 


0l0x5 




































05X10 


O5XIO 


c 4 


O5X5 




05x10 


O5X10 


05x5 


05x5 


\O5xlO 


O5XIO 


05x5 


05x5 / 




\ 05x10 


O5X10 


05x5 


C 4 y 



(A.l) 



19 



w = (l,2,2,2,l,2,2,l,2,l,0 lx20 ) T , 
y 3 = (2, 0, 2, 0, 2, 0, 0, 0, 0, 2, 0, 2, 2, 0, 0, 0, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) T , 
y A = (0, 2, 2, 0, 0, 0, 2, 2, 0, 0, 0, 0, 0, 0, 0, 2, 2, 0, 2, 2, 2, 0, 2, 0, 2, 0, 0, 0, 0, 0) T , 
y 5 = (2, 2, 0, 2, 2, 2, 0, 2, 0, 2, 0, 0, 0, 0, 0, 4, 2, 0, 2, 0, 2, 2, 0, 2, 0, 0, 2, 0, 0, 0) T , 
y 6 = (4, 2, 0, 2, 0, 2, 2, 0, 2, 0, 0, 2, 0, 0, 0, 4, 2, 2, 2, 4, 2, 0, 0, 2, 2, 0, 0, 0, 0, 0) T , 
y 7 = (4, 2, 2, 2, 4, 2, 0, 0, 2, 2, 0, 0, 0, 0, 0, 4, 4, 4, 2, 0, 2, 2, 0, 2, 0, 0, 0, 0, 0, 2) T , 
y 8 = (4, 4, 4, 2, 0, 2, 2, 0, 2, 0, 0, 0, 0, 0, 2, 6, 4, 4, 2, 4, 2, 0, 4, 2, 2, 0, 0, 0, 0, 0) T , 
y 9 = (6, 4, 4, 2, 4, 2, 0, 4, 2, 2, 0, 0, 0, 0, 0, 8, 4, 4, 2, 0, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0) T , 
y 10 = (8, 4, 4, 2, 0, 4, 4, 4, 0, 0, 0, 0, 0, 0, 0, 8, 4, 6, 4, 8, 2, 0, 4, 2, 4, 0, 0, 0, 0, 0) T , 
y n = (8, 4, 6, 4, 8, 2, 0, 4, 2, 4, 0, 0, 0, 0, 0, 8, 6, 6, 2, 0, 2, 6, 4, 2, 0, 2, 0, 2, 2, 0) T , 
y 12 = (8, 6, 6, 2, 0, 2, 6, 4, 2, 0, 2, 0, 2, 2, 0, 10, 6, 4, 4, 8, 2, 0, 4, 2, 4, 0, 0, 0, 0, 0) T , 
2/i3 = (10, 6, 4, 4, 8, 2, 0, 4, 2, 4, 0, 0, 0, 0, 0, 12, 6, 4, 4, 0, 6, 6, 4, 2, 0, 0, 0, 0, 0, 0) T , 
y u = (12, 6, 4, 4, 0, 6, 6, 4, 2, 0, 0, 0, 0, 0, 0, 10, 6, 8, 6, 12, 4, 0, 0, 4, 4, 0, 0, 0, 0, 0) T , 
y 15 = (10, 6, 8, 6, 12, 4, 0, 0, 4, 4, 0, 0, 0, 0, 0, 8, 10, 6, 6, 0, 4, 8, 4, 4, 0, 2, 2, 0, 0, 0) T , 
and introducing the recurrence relation 

V2n = F Q y n , y 2n +i = F x y n , n>8 
one has the relation [8]: 

u n+1 = w T y n . (A. 2) 

We finally introduce two new matrices in ]R 20x20 that rule the asymptotics of 
u n : 

A = 10X10 ,A 1= 1 . (A.3) 

\Di B 1 j V0 10 xioC ly / 



B Proof of Theorem 2 



In this appendix we give a proof of Theorem 2. We first find a common in- 
variant cone K for the matrices A ,Ai. We then show that the products 
F{n) = ■ ■ ■ Fd k are asymptotically equivalent to their corresponding prod- 
uct A(n) = Adj_ ■ ■ ■ Ad k . We then shed some light on the vectors z n , the restric- 
tion of y n to M 20 : their norms can be considered as norms of the products 
A(n). We finally show that WA^ ■ ■ • Ad h \\ is equivalent to \\Adj_ ■ ■ ■ Ad k _ 4 \\. We 
end with the proof of Theorem 2 that puts all this together. 



20 



Let us first establish a special property of the matrices Ao,A±. Consider the 
following sets: 

P= (x G R 20 ,x > 0, Xi > for all i <£ {5,10,17}}, 



Q = Ix G R 20 ,x> 0, Xi > for all % £ (7, 15, 20} V 

Let S = PUQ. For any e > let p £ G R 20 , (p E )i = -e for i G {5, 10, 17} and 
{p e )i = 1 otherwise, let also q £ G M 20 , {q e )i = —£ for f G {7, 15, 20} and (q E )i = 
1 otherwise. This is easy to verify by direct calculation that there is e > 
such that A p £ , A q e G P and A\q E , A 1 p £ G Q (this is the case for instance for 
e = 1/4). Take any such e and denote K = {v + tp £ + sq £ , v G M 2 ^, t, s > 0}. 
Clearly, K is a convex closed pointed cone. 

Lemma 2 VFe have AiK C if and AjiS 1 C S* /or i = 0, 1. Moreover, S C intif. 

Proof. Let x E K : x = v + tp £ + sg £ , t> > 0, t, s > 0. Since A;t> > 0, A,ip £ > 
and Ajg e > we see that AiX > 0. Thus, AiK C Ji . 

Further, for an arbitrary i e S let x m i n = minxj. Assume x G P (the 

proof for the case x G Q is literally the same). Since x > x m ; n po w e have 
Aqx > x min A p > x min p G P and Aix > x min Aip > a; min g G Q. Therefore 
AiX G S. 

Finally, for any x G S and y G M 20 such that \y\oo < min{|£x min , |e} we 
have x + y E K. This proves that x G intif. □ 

Corollary 5 For any x G S and for any sequence d\, . . . ,dk we have A^ ■ ■ ■ Ad k x G 
intK. 

Lemma 3 Suppose n > 1 and consider the binary expansion c4 . . .d\ of the 
number n — 1. We define A{n) = A^ ■ ■ ■ Ad k and similarly for F(n). Then, 
for any matrix norm one has: 

||F(n)||x|L4(n)||. 

Proof. Since all matrix norms are equivalent, it suffices to consider any norm. 
Obviously \\F dl ■ ■ ■ F dk \\i > \\A dl ■ ■ ■ A dk \\i (because the matrices Ai are sub- 
matrices of Fi), hence it remains to prove the opposite inequality: there is a 
positive constant C such that \\A dl ■ ■ ->-Lj fc ||i > C||F dl • • --F^JIi for all k G N 
and d\, . . . ,dk. We consider the case where d^ = 1; the proof for the other 
case is similar. Let m < k — 1 be the biggest number such that d m = 0. If the 
sequence has no zero, we fix m = 0. Let Ai, Hi, and Ri denote respectively 
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the upper left, the upper right and the lower right corners of the matrix Fi in 
bloc representation A.l. Then the product Fd 1 . . . F dk has the following form: 
the left upper block is A dl ■ ■ ■ Ad k , the left lower block is zero, the right lower 
block is Rd 1 ■ ■ ■ Rd k , finally, the right upper block is 



p =i s=i 



n ^ 

j=p+i 



(B.l) 



By convention, the product over an empty set is one. Since RqR± = 0, the 



right lower block is zero and block (B.l) becomes (lTj=i Ai, ) Hd p R 1 

p=m ^ ' 

whose norm can be estimated from above as 



k—p 



k p-1 
p=m j=i 



k—p 



(B.2) 



where H = max{||i?o||i, ||-^i||i}- It was shown in [8] that the sum of entries 
of the matrix R\ does not exceed C2 l for any i = 0, 1 and I > 1, where C > 
is a constant. Hence ||-Ri _p ||i < C2 k ~ p . Thus, 



< 



.4 



EC Y, 2fc " P II A * 



p-1 



fB.3) 



On the other hand, for any p > m we have \A^ ■ ■ ■ Ad k l\ = \A^ 1 ■ ■ ■ Ad p _ 1 {Ad p A\~ 
where 1 is the vector of ones. By Corollary 5 the vectors Aj AJ1 belong to 
int-fT for all r G N. Moreover, the vector Ajl/|-i4.£l| converges to v± (the 
Perron- Frobenius eigenvector of A±) as r — > oo and V\ G S, hence V\ G iritK 
(Lemma 2). Therefore there is a constant C\ > such that ^(A^A^l) > G\ 
for all r G N. Applying now Lemma 1 for x — A d Ai~ p l, we get 



1-4 



di ■ ■ • A dk I 



> 



C 2 \A dl 



A dp _ 1 {A dp A\~ p l) 



> CadHAt ■■■A dp _ 1 \\ ■ \A dp A\- p l\ > CsX^WAd^'-Ad^l 

where A = p(Ai) (by the same reasoning we have \Ad A[l\ > 7(^1) \\Ad \ \ ■ 
\A{1\ > CX'\ indeed, A\\ G intK, and thus |L4 d *v^l|| x \\A d ?\\ ■ \A\\\). 
Thus, ||A dl ■ ■ ■ Aip-ill < C^\ p ~ k \A dl ■ ■ f o r all p < k. Substituting this 

in (B.3) and taking into account that < 1 (because A = p(A) > 2.42) 
we take the sum of the geometrical progression and get ■ • ■ F dk ||i < 
C4IIA21 ' ' ' AzJ|i, where C4 is some constant. This concludes the proof. □ 



Lemma 4 For any n we have u n+ i < 2u n . 

Proof. If a word of length n + 1 is overlap-free then so is its prefix of length 
n. On the other hand, at most two overlap-free words of length n+ 1 have the 



22 



same prefix of length n. 



□ 



Lemma 5 Let the vectors y m G M 30 be the solution of the recurrence equation 
4, and z m G M 20 be the vector with the first 20 entries of y n . We have z m G S 
for each m = 64, . . . , 127. 

Proof. The proof is by direct calculation. □ 



Lemma 6 Suppose n G N and df- . . . d\ is the binary expansion of n — 1; then 
\\A{n)\\ x \\A\n)\\, where A(n) = A dl ■ ■ ■ A dk and A'{n) = A dl ■ ■ ■ A dk _ 4 . 



Proof. The inequality 



A dl ■ ■ ■ A a 



tiplicativity of the norm. For the other direction, we have: 



< C 



A d ,--- A dh 



is obvious by submul- 



A dl ■ ■ ■ A a 



Ah,--- A/,. 1 



A dl ■ ■ ■ A dk _ 4 (A ( 



: fe _ 3 • • • A dk l / 



(B.4) 



Corollary 5 yields A dh _ 3 . . . A dk l G 'miK for all dks, ■ ■ ■ , d^ G {0, 1}. Applying 
now Lemma 1 we get 



6 = * m j n c /n i x 7(A dk _ 3 ...A dk l)> 0. 
<4_ 3 ,...A€{o,i} 



Therefore, for some C\ > 0, 



A dl ■ ■ ■ A dk _ 4 



< S 
< 



A dl ■ ■ ■ A dk _ 4 



A dk _ 3 A dk _ 1 A dk l 



A dl ---A d ,l 



d k - 



(B.5) 



Combining this with (B.4) we get 
We are now able to prove Theorem 2: 



A dl ■ ■ ■ A dk _ 4 



< a 



M---A dk 



□ 



Proof of Theorem 2. Let g be the vector of M 30 , whose first 20 entries are ones 
and the last 10 entries are zeros. Let also m = c4_ 3 + c4-22 + 2 2 c4_i + 2 3 dk- 
Since w < 2g, we have 



u n < 2(y n _ u g 
< C \\F dl -.-F t 



2[F dl - ■■F dk 4 y m ,g 



\A dl ■ • ■ A c 



(B.6) 



where Cq does not depend on n (the first two relations are direct from funda- 
mental assertions A. 2, the third relation comes from the fact that y m and g 
are bounded, and the last equivalence is by Lemma 3). Combining Lemma 6 



and (B.6) gives u n < Cj 



A 



di ■ ■ 



■A dk 
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Let us now prove the opposite inequality. Lemma 4, together with the fact that 
by construction, the first ten entries of y n are equal to the entries 11, . . . , 20 
of y n _i, imply that u n > ~(u n + u n+ xj > \{y n ^gy Furthermore, for n > 2 7 

we have (y n -u g) = [F dl ■ ■ ■ F dk7 y h , where I = £j =0 4-e+j2 J . Thus, 



u, 



> -{F dl ---F dk _ 7 y h g 



(B.7) 



On the other hand, defining z\ G M 20 as the vector with the first 20 en- 
tries of y t : (F dl ■ ■ ■ F dk _ 7 y h gj > (A dl ■ • • A dk _ 7 z u l)= A dl ■ ■ ■ A dk _ r z t j . By 
Lemma 5 we have z\ G mtK for all I G {64, . . . , 127}, and we can define 



■A 



d k _ 7 Zl 



h = min 7(2/) > such that 

64<Z<127 V ' 

Ca = min \zi\i. Combining this with (B.7), we obtain 

64<K127 



>hC 4 \\A dl ---A a 



u n > C 5 \\A dl ■ ■ ■ A dk _ 7 \ 
Now, by submultiplicativity of the norm, 

U n > C 6 \\A dl •••A d J| 



where 



(B.8) 



□ 



C The ellipsoidal norm 



Define 



Pi = 



31.3 


7.5 


2.3 3.3 


7.5 


57.7 


10 6.3 


2.3 


10 


59.9 11.3 


3.3 


6.3 


11.3 48.5 


-0.4 


18.4 


0.4 4.6 


-0.3 


35 


29.2 13.5 


0.3 


16.3 


4.2 10.8 


0.4 


-5.8 


10.1 2 


3.7 


13.8 


8.2 6.9 


0.3 


5 


0.8 1 



-0.4 -0.3 0.3 0.4 3.7 0.3 

18.4 35 16.3 -5.8 13.8 5 

0.4 29.2 4.2 10.1 8.2 0.8 

4.6 13.5 10.8 2 6.9 1 

36.4 23.5 22.6 4.4 8.9 -1.2 

23.5 105.9 38.4 9.5 33.7 6.1 

22.6 38.4 59 2.7 17.4 9.2 
4.4 9.5 2.7 38.6 14.8 -1.7 
8.9 33.7 17.4 14.8 57.5 8.6 

-1.2 6.1 9.2 -1.7 8.6 42.3 y 
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p 2 = 



— 10.4 


— 1.7 


— 18.1 


—0.4 


—5.8 


—5.1 


—4.9 


—0.8 


-2.7 


—0.9 


— 11.1 


—22.4 


—8.2 


— 14.7 


—9.9 


—30.3 


— 16.7 


— 11.3 


— 16.9 


—6.6 


— 2.2 


—16.4 


— 15.8 


—5 


—8.5 


—7.2 


r* A 

—5.4 


— 18.5 


—3.5 


—3.4 


—0.2 


—13.6 


—5.2 


—9 


— 10.7 


— 14.6 


—9.2 


— 1.6 


— 11.3 


— 1.1 


A f 

—4.6 


—17 


— 13 


—9.1 


—0.6 


— 11.2 


— 23.9 


—7 


— 12.1 


0.3 


—5.9 


—26.4 


—27 .4 


— 17.4 


—31 


—37.6 


—28 


—4.4 


— 27.3 


T A 

—7.4 


-1.4 


-19.3 


-11.6 


-10.8 


-22.3 


-17.9 


-11.7 


-11.3 


-12 


-9.8 


-6.3 


2.1 


1.7 


-3.4 


3.2 


-7.6 


0.2 


-5.2 


-3.1 


-1.4 


-7.4 


-15.9 


-4.7 


-6.7 


-12.2 


-17.3 


-11.6 


-5.3 


-6.8 


-1.6 


1.3 


-5.7 


-3.6 


-3.2 


-0.4 


-6.1 


-9 


-1.4 


-6.9 


0.4 



p 4 = 



29.1 8.3 


-1.6 


4.8 


-1.3 


-4.4 


0.6 


1.7 


7.5 


1.1 


8.3 47.3 


13.6 


2.8 


11.7 


19.8 


17.4 


0.6 


10 


3.7 


-1.6 13.6 


46.6 


10.4 


6.5 


24.9 


11.8 


6.5 


12.5 


1.4 


4.8 2.8 


10.4 47.6 


5.1 


8 


7.6 


5.1 


3.7 


1.8 


-1.3 11.7 


6.5 


5.1 


32.8 


19.5 


19.4 


7.6 


6.7 


-0.2 


-4.4 19.8 


24.9 


8 


19.5 


64.8 


16.2 


11.4 


13.8 


6.8 


0.6 17.4 


11.8 


7.6 


19.4 


16.2 


56.7 


7.6 


12.2 


6.5 


1.7 0.6 


6.5 


5.1 


7.6 


11.4 


7.6 


38.7 


11.2 


-1 


7.5 10 


12.5 


3.7 


6.7 


13.8 


12.2 


11.2 


55.6 


4.2 


1.1 3.7 


1.4 


1.8 


-0.2 


6.8 


6.5 


-1 


4.2 


43.8 



p=M. 

Then one has the relations: 

A 1 PA - (2.5186) 28 P -< 0, VA e S 14 . 
As explained in [3] this suffices to prove that p(S) < 2.5186. 
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D The vector x 

Define 

x = (153, 0, 60, 0, 50, 56, 99, 0, 58, 1, 157, 81, 0, 113, 0, 72, 0, 99, 0, 0) T . 
Then one has the relation 

B(Ax - rx) > 0, WB e S 6 , A e S 16 , 
x > 0, 

with r = 2.41 16 . this proves that > 2.41. 
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